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Measurements by the Belle Collaboration of the exclusive production of two charmonia in e + e~ 
annihilation differ substantially from theoretical predictions. Till now, no conclusive explanation for 
this remarkable discrepancy has been provided. Even the origin of the discrepancy is not identified, 
yet. We suggest that the measurement of four-charm events in Belle data must provide a strong 
constraint in identifying the origin of this large discrepancy. Our prediction of the cross section for 
e + e~ — > cccc, in lowest order in strong coupling constant, at ^/~s = 10.6 GeV is about 0.1 pb. If 
measured four-charm cross section is compatible with the prediction based on perturbative QCD, it 
is very likely that factorization of hadronization process from perturbative part may be significantly 
violated or there exists a new production mechanism. If the cross section for the four-charm event 
is also larger than the prediction like that for the exclusive J/ip-\-rj c production, perturbative QCD 
expansion itself will be proved to be unreliable and loses predictive power. 

PACS numbers: 13.66.Bc, 12.38.Bx, 14.40.Gx 



The Belle Collaboration has measured the cross sec- 
tion for J/iJj + rj c by observing a peak in the momentum 
spectrum of the J/ip that corresponds to the 2-body final 
state J/ip + T] c 1] . The measured cross section is 



a[J/iP + i lc ] x B*°[> 4] = 



± 9) 



fix 
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where B 71 " [> 4] is the branching fraction for the r\ c to de- 
cay into at least 4 charged particles. Since B Vc [> 4] < 1, 
the right side of Eq. |Q is a lower bound on the cross 
section to produce J/ip + rj c . This lower bound is about 
an order of magnitude larger than the predictions by 
Braaten and Lee [2|, and by Liu, He, and Chao 3j of non- 
relativistic QCD(NRQCD) [4| in the nonrelativistic limit. 
The cross section was calculated previously by Brodsky 
and Ji |5| using perturbative-QCD factorization formal- 
ism 0. But they did not give an analytic expression for 
the cross section. Recently, Brodsky, Ji, and Lee redid Q 
the calculation given in Ref . |f| ■ They found exact agree- 
ment Q with the result based on NRQCD HQ. Cur- 
rently, the cross section for the process e + e~ — » J/ip + r] c 
shows the largest discrepancy between theory and data 
available within standard model. This presents a chal- 
lenge to our current understanding of charmonium pro- 
duction based on perturbative QCD framework. 

A few theoretical studies have been carried out in or- 
der to explain the large discrepancy. Bodwin, Braaten, 
and Lee proposed 0, that the Belle data for J/ip + rj c 
may include the J/ip + J/ip events because the width of 
the r\ c peak in the recoil mass distribution for inclusive 
J/ip production measured by the Belle Collaboration is 
wide enough to accommodate J/ip events. The large en- 
hancement from photon fragmentation in the two-photon 



mediated process e + e~ — > J/ip + J/ip overcomes the 
suppression factor a 2 /cv 2 in couplings compared to the 
e + e~ — ► J/ip + T] c [alPj- Brodsky, Goldhaber, and Lee 
introduced an exotic scenario that the Belle J/ip + r\ c 
signal may include the associate production of J/ip and 
spin- J glueball Qj, J = 0, 2 [Ic}. The Belle Collaboration 
carried out an updated analysis 01 motivated by these 
proposals. According to the Belle analysis, no events 
for J/ip + J/ip have been detected and the upper limit 
for cr(e+e~ -> J/ip + J/ip) x B J ^[> 2] is 9.1 fb at the 
90% C.L., which is consistent with the prediction given in 
Refs. H,@- Here B J ^[> 2] is the branching fraction for 
the J / ip to decay into more than 2 charged particles. The 
measured cross section for e + e~ — > J/ip + J/ip, however, 
does not explain the large fraction of the J/ip + rjc sig- 
nals. The Belle Collaboration also analyzed the angular 
distribution of J / ip in order to identify if the data include 
the associated J/ip + Go signals. The predicted angular 
distributions are proportional to cos 2 9 for J/ip + rj c 
and J/ip + G 2 E3. and sin 2 for J/ip + G 0, where 6 is 
the scattering angle of the J/ip in the e + e~ cm. frame. 
The updated analysis show the measured distribution is 
proportional to cos 2 0, which ruled out the spin-0 glueball 
scenario 0] . Since the angular distribution for J/ip + Q2 
has the same form as that for J/ip + i] c , spin-2 glueball 
has not been ruled out, yet. 

Another scenario is that higher-order corrections in 
strong coupling a s may be huge [l^. If it is true, pertur- 
bative expansion is not a proper method to predict the 
cross section. If it is not, it is probable that the factor- 
ization of long-distance factor involving hadronization is 
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FIG. 1: Two of the 8 topologically distinct Feynman 

seriously violated or there exists a unknown production 
mechanism which we do not understand, yet. Once we 
can estimate the size of the perturbative QCD correc- 
tions to this process, it might be easier for us to iden- 
tify the origin of this large discrepancy. Unfortunately, 
the next-to-leading-order corrections to the cross section 
of exclusive J/ip + r\ c process is not available. A com- 
prehensive review on recent developments in quarkonium 
physics can be found in Ref. [l3|. 

In this paper, we introduce an economical method to 
check if the perturbative QCD corrections to the cross 
section is indeed large enough to explain the discrep- 
ancy. If we consider inclusive four-charm-hadron produc- 
tion, the prediction for the cross section can be expressed 
as inclusive cccc production rate <r(e + e~ — * cccc + X). 
This is analogous to estimating a(e + e~ — ► hadrons) 
by J2 q a ( e+e ~~ ~~ * <Z<Z + We expect er(e + e~ — > 

cccc + X) « cr(e + e~ — » cccc) is a good approximation 
at sfs — 10.6 GeV. Unlike the prediction for J/tp + r\ c 
cross section, the prediction for the inclusive four-charm- 
hadron production rate purely consists of short-distance 
factor. Corresponding long-distance factor for hadroniza- 
tion is of order 1. Since this process involves the same 
Feynman diagrams for exclusive J/ip + Tj c production, 
the measurement of the cross section for four-charm- 
hadron production will provide an important informa- 
tion in estimating the size of the short-distance coeffi- 
cient for J/ip + r) c cross section. We present our predic- 
tion for inclusive four-charm-hadron production by cal- 
culating cr(e + e~ — > cccc + X) in order a 2 a 2 , which is 
leading order in strong coupling constant. If our leading- 
order prediction is comparable to the measured value, it 
is very probable that the QCD higher-order corrections 
to the J/tp + i] c cross section is small. Then the large 
discrepancy in J/ip + rj c cross section may be due to the 
violation of factorization or existence of new production 
mechanism. If the measured cross section for the four- 
charm-hadron inclusive production is much larger than 
our prediction like the case of J/ip + ry c , it is very likely 
that perturbative QCD corrections to J/ip + r/ c cross sec- 
tion is large enough to explain the discrepancy, which 
leads to the failure of reliability in perturbative expan- 
sion. 

In leading order in strong coupling a s , cccc can be pro- 
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duced at order a 2 a 2 . There are two topologically distinct 
Feynman diagrams generating two pairs of cc, which are 
shown as Ali and M 2 in Fig- Wi & ) an d H^b), respec- 
tively. Momenta for the involving particles are assigned 
as e~ (fci)e + (A:2) — > c{pi)c(p 2 )c(p 3 )c(p4). The amplitude 
for the two diagrams shown in Fig.n& re 

. . (47r) 2 e c aa s _ 

-iMi = I—. ; r^V e (k 2 )j a U e (ki) 

S(P2 +P3) 

x u{p 3 )T a i P vfa)u{p l )T a Hfv{p A ), (2) 

where s — (ki + k 2 ) 2 , e c = | is the fractional electric 
charge of the charm quark, and a is the SU(3) color index 
for the virtual gluon. The vector indices a and /3 are for 
the virtual photon and gluon, respectively. We suppress 
the spin and color indices of the charm quarks in Eq. J5J. 
For i = 1 or 2 the tensors H" 13 in Eq. J2J), which are 
matrices in spinor space, are defined by 

Hf = ^A( Pl +p 2 +p 3 ) 1 a , (3a) 
Hf = 7 Q A(-p 2 -p 3 -p 4 )7 /3 , (3b) 

where A(p) = (]/> + m c )/(p 2 — m 2 ). 

There are 6 more Feynman diagrams that can be ob- 
tained from the two amplitudes Mi and M 2 by exchang- 
ing two charm quarks and two antiquarks, respectively, 
as 

M 3 = -JWMi, M 4 = -Pi^ 3 M 2 , 

M 5 = -P 2 ^4Mi, M G = -P 2 ^M 2 , (4) 

M 7 = +Px^zP 2 ^aMi, M s = +Pi^ 3 P2^iM 2 . 

where Pi^j is the operator exchanging two particles with 
momentum indices pi and pj shown in Fig. ^ The signs 
of M 3 through Ms in Eq. 10} are determined by the an- 
tisymmetricity of Fermi statistics in exchanging identical 
fermions among the final-state particles. 

The total cross section for the process is expressed as 

where M = T* =1 Mi, and the factor (2!) 2 in Eq. © 
is divided in order to avoid double-counting of identi- 
cal final-state particles. The summation notation 
in Eq. © stands for averaging over initial spin states 
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FIG. 2: Total cross section <r(e + e — > cccc) at ^/s = 
10.6 GeV in pb as a function of m c , where a — 1/137 and 
a s = 0.2. 



and summation over final color and spin states. The 
four-body phase-space element g?<I>4 in Eq. JSJ can be 
parametrized by 
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where my is the invariant mass of pi + pj and dQ*^ is the 
solid angle element of pi + pj in the rest frame of pi + pj . 
Their physical regions are 2m c < mi3 < y/s — 2m c and 
2m c < m24 < y/s— mi3. The three- momenta pi ofpi and 
P2 of P2 are defined in p\ + p% and P2 + Pa rest frames, 
respectively. The three-momentum P and solid angle 
element dVl are for p\ + P2 in the e + e~ cm. frame. In- 
tegrating the differential cross section (0 over the phase 
space ©, we get the total cross section for e + e~ — » cccc. 
We compute the J2\M\ 2 in Eq. © using REDUCE 
and carry out the phase-space integral in Eq. © making 
use of the adaptive Monte Carlo routine VEGAS 0] . As 
a check, we carry out the same calculation using Com- 
pHEP 16] . Our analytic result for X)l-^| 2 an d numer- 
ical values for the total cross section agree with those 
obtained by using CompHEP. 

Our predictions for the inclusive four-charm-hadron 
cross sections in e + e~ annihilation at ^fs =10.6 GeV de- 
pending on the charm-quark mass m c is shown in Fig. [21 
The cross section for e + e _ — > cccc is very sensitive to the 
value of m c . For a = 1/137, a s = 0.2, m c =1.5 GeV 
cr(e + e _ — > cccc) =97 fb. The cross section varies from 
0.31 pb at m c =1.2 GeV to 24 fb at m c =1.8 GeV. The 
cross section decreases as m c increases mainly because 
available phase space shrinks. If one can increase the cm. 
energy of the e + e~, the m c dependence will decrease. In 
previous analyses for exclusive two-charmonium produc- 
tion cross sections, the next-to-leading order pole_mass 
m c = 1.4 ± 0.2 GeV has been used for the m c 



FIG. 3: Differential cross section da/dm cc in fb/GeV with 
respect to the invariant mass m cc = mra of cc for e + e~ an- 
nihilation into cccc, where m c — 1.5 GeV, a = 1/137, and 
a s — 0.2. Physical range of the m cc is from 2m c to y / s-2m c . 
The area under the curve is the integrated cross section 97 fb. 



However, the cross section for e + e~ 
sitive to the charm-quark mass m c . 
cross section of order a 2 is cr(e + e~ - 
relative errors of 3 x 10~ 3 for m c = 



— > cc is not sen- 
The lowest-order 
cc) = 1.0 nb with 
1.5 ± 0.3 GeV at 



= 10.6 GeV. In Ref. [Fjj the total cross section for 
e+e~ — > cccc at y/s = 10.6 GeV is predicted. If we use 
the input parameters a s = 0.24 and m c = 1.4 GeV given 
111 Ref. hi, we get 0.210 pb, which is different from the 
prediction 0.237 pb given in Ref. 17] by about 13%. 

In Fig. 13 we show the differential cross section with re- 
spect to the invariant mass of cc. This is the prediction 
for da(e + e~ — > cc + X)/dm cc in leading order in a s . Ex- 
perimentally, this differential cross section can be com- 
pared with the J2 H H , da{e + e~ — > HH' + X)/drriHH', 
where H and H' are charm hadrons, which do not in- 
clude anticharm. 

Finally, we estimate the number of four-charm-hadron 
events that could be detected by the Belle Collabora- 
tion. Production rate for baryonic states such as A c will 
be small and we do not include the contribution in the 
following rough estimate. Based on heavy-quark spin 
symmetry, the relative rates for a c quark fragmenting 
into the charm mesons are D + : D° : D + * : D°* = 
1:1:3:3. Because the two spin-triplet states de- 
cay into spin-singlet states by 100% with branching frac- 
tions Br[D+* -> J D°7r+]=70%, Br[£>+* -> £>+7r°]=30%, 
Br[L>°* -► L> 7T°]=62%, and Br[L>°* -► L>°7]=38%, we 
may only consider charm meson pairs made of either D + 
or D°. Resulting fragmentation probabilities are approx- 
imately P[c -► D+ + X] pb \ and P[c -> D° + X] ps |, 
respectively. Therefore, cr[e + e _ — * D + D + + X] s» 
* cc + X], a[e + e- D+D° + X] ps 

cc + X], and a[e+e- -> D°D° + X] w 



4(r]e+e- 



16 



4 



cc + X]. The detection rate will suffer losses 
from branching fractions Br[D + — > K ~n + ir + ] = 9.2% 
and Br[_D° — > ir~7r + ] = 3.8%, and detection accep- 
tance/efficiency ss 80% for each charged particle in the 
decay products of D + or D° . With er[e + e~ — > cc + X] w 
0.1 pb and current integrated luminosity £ sa 300 fb _1 
we expect roughly 30 events will be detected by the Belle 
detector. Even if we consider the uncertainties from a a 
and m c in our prediction, we expect at least about 10 
events will be detected by the Belle Collaboration. If 
there is a large QCD corrections, the number of events 
will be increased into several hundreds. 

In summary, we have calculated the cross section for 
e + e~ annihilation into cccc. Assuming quark-hadron 
duality, the cross section for the inclusive four charm 
hadrons is predicted to be about 0.1 pb. The com- 
parison of this prediction with the measured cross sec- 
tion for the four charm hadrons at B-factories will pro- 
vide a strong constraint in determining the origin of the 
large discrepancy between prediction and Belle data for 
exclusive J/ip + rj c production in e + e~ annihilation at 
y/s = 10.6 GeV. The measurement will also provide a 
useful information in explaining large cross section for 
J/ip + cc + X measured by the Belle Collaboration jlj 
compared to the NRQCD predictions (l8j . 
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